clf;
data = load('gp_data');
fns = load('gp_function');
n_class =  3;
n = size(fns,1)/n_class;
fns = reshape(fns, n, n_class);
hold on;
plot(1:n, fns(:,1), '-1;f 1;');
hold on;
plot(1:n, fns(:,2), '-2;f 2;');
hold on;
plot(1:n, fns(:,3), '-3;f 3;');

% yes I know line graphs are inappropriate here
% they're still easier to read
% or at least they are with no errorbars

% try a variant of
% errorbar (x, f(x), ym, yp);
% to get errors. More helpful-like

% NTS: colours are indexed rgbmcw